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Broadside  Radar  Cross  Section  of  the 
Perfectly  Conducting  Cube 


1.  INTRODUCTION 

With  its  six  faces,  twelve  edges,  and  eight  corners,  the  cube  presents  a 
challenging  problem  for  the  theories  and  measurement  of  electromagnetic  scatter¬ 
ing.  The  magnetic  field  integral  equation  (MFIE)  solution  of  Tsai,  Dudley  and  Wilton, 

for  backscattering  cross  section  from  a  perfectly  electrically  conducting  (PEC)  cube 

2 

under  broadside  plane-wave  incidence,  agreed  well  with  the  measurements  of  Ryan 
up  to  a  cube  side  length  (s)  of  about  0.  5  wavelengths  (X).  For  s  beyond  about  0.  5X 
the  spurious  solutions  admitted  by  the  MFIE  begin  to  contaminate  the  numerical 
solution,  and  comparisons  between  theory  and  experimental  measurements  become 

3 

unreliable.  To  overcome  this  spurious  resonance  problem,  Yaghjian  introduced 
"augmented”  integral  equations  and  applied  the  augmented  magnetic-field  integral 
equation  (AMFIE)  to  the  problem  of  broadside  backscattering  from  the  conducting 
cube.  Because  his  numerical  solution  did  not  exploit  the  symmetry  of  the  cube, 

(Received  for  publication  26  January  1984) 

1.  Tsai,  E.L.,  Dudley,  D.G.,  and  Wilton,  D.  R.  ( 1 974)  Electromagnetic  scatter¬ 

ing  by  a  three-dimensional  conducting  box,  J.  Appl.  Phys.  ,^(10:4393-4400. 

2.  Ryan,  ('.  E. ,  Jr.  (1970)  Diffraction  Analysis  of  Scattering  by  a  Cube  With 

Application  to  the  Time  Response  Waveforms,  The  Ohio  State  University 
Electroscience  Laboratory  Report  2415-3. 

3.  Yaghjian,  A.  D.  ( 198  1 )  Augmented  electric-  and  magnetic -field  integra l 

equations.  Radio  Sci.  ,  16(6):9  87  - 100 1.  (See  also  RADC  In-house  Report 

TR  -81  -45  under  the  same  title,  A  DA  103946.  ) 


and  his  matrix-inversion  subroutine  required  that  the  matrix  be  stored  on-line, 
computer  storage  capacity  limited  the  accurate  solution  of  the  AMFIE  to  a  cube  side 

length  of  about  1A. 

2  4 

Ryan  '  applied  the  geometrical  theory  of  diffraction  (GTD)  to  determine  the 
high  frequency  broadside  backscattered  fields  from  conducting  cubes  at  harmonically 
related  frequencies  which  were  combined  with  low  frequency  experimental  data  and 
transformed  to  obtain  time-domain  responses.  Although  no  radar  cross  section 
(RCS)  vs  frequency  curves  were  shown  in  Reference  2  or  4,  a  similar  high  fre¬ 
quency  diffraction  (HFD)  solution  for  broadside  RCS  is  extracted  from  the  enhanced 

3|C 

high  frequency  solution  of  Section  2.7.  Comparison  between  the  HFD  curve 
(Figure  3)  and  the  AMFIE  curve  of  Reference  3  reveals  that  there  is  a  poor  agree¬ 
ment  between  the  HFD  and  AMFIE  solutions  within  the  region  of  accuracy  (s  less 
than  about  (<)  A]  of  the  AMFIE  solution.  Also,  as  discussed  in  Section  2.  7, 
inclusion  of  higher  order  edge  diffracted  rays  in  the  HFD  solution  did  not  improve 
the  agreement  in  this  region  nor  enhance  the  resonance.  Moreover,  for  s  >  A  the 
HFD  backscattering  solution  differs  slightly  from  the  physical  optics  (PO)  solution 
which  is  plotted  in  Figure  4.  Thus  we  have  little  assurance  that  the  conventional 
HFD  or  GTD  solution  is  a  significant  improvement  over  the  physical  optics  solution 
for  broadside  RCS  vs  frequency  of  a  cube. 

The  purpose  of  the  present  work  is  to  predict  theoretically  the  broadside  back- 
scattering  cross  sections  of  perfectly  conducting  cubes  with  reasonable  engineering 
accuracy  (*  ±  1  dB)  from  arbitrarily  low  to  arbitrarily  high  frequencies,  and  to 
compare  the  results  to  measured  RCS  data  covering  a  frequency  range  that  extends 
well  beyond  the  resonance  region.  We  take  a  composite  approach  starting  with  an 
MFIE  solution  similar  to  that  of  Reference  1  for  s  less  than  0.  4  A.  For  s  greater 
than  0.  4  A  the  AMFIE  is  solved  utilizing  x-y  symmetry  in  order  to  reduce  the  number 
of  unknowns  and  equations  by  a  factor  of  four,  thereby  allowing  an  accurate  numerical 
solution  on  the  available  Cyber  7  50  computer  up  to  a  side  length  of  about  1.5  A. 

Beyond  1.5  A  an  "enhanced  high  frequency  diffraction"  (EHFD)  solution  is  developed 
that  integrates  the  currents  induced  on  the  leading  face  and  edges  of  the  cube  to 
determine  the  fields  impinging  upon  the  trailing  edges.  The  mutual  validity  of  the 
EHFD  and  integral  equation  solutions  is  ascertained  by  their  close  agreement  in  the 
intermediate  frequency  range  0.  4  A  1  s  <  1.5A.  (Also,  the  EHFD  solution  reveals 

4.  Ryan,  C.  E. ,  Jr.  (1971)  Time- response  waveforms  for  a  cube  from  measured 
data  and  diffraction  analysis.  Radio  Sci.  6(8,  9):801-804. 

“  - -  WN 

The  HFD  solution  differs  from  the  GTD  solution  of  References  2  and  4  because  the 
HFD  solution  includes  doubly  diffracted  ravs  (using  the  term  "doubly  diffracted" 
as  defined  in  Section  2.  5),  and  because  in  the  limiting  procedure  used  to  obtain 
broadside  RCS  from  the  GTD  diffraction  coefficients,  the  fringe  current  fields 
fron  the  front  face  are  not  entirely  recovered. 


clearly  the  dominant  resonance  and  shows  decent  agreement  with  the  integral  equation 
solution  throughout  the  resonance  region,  s  <  0.  5  X  . )  The  composite  RCS  vs  fre¬ 
quency  curve  agreed  to  within  about  ±  1  dB  with  measured  data  taken  at  X-band  with 
aluminum  cubes  ranging  in  side  length  from  0.  15  to  4  wavelengths. 


2.  THE  ENHANCED  HIGH  FREQUENCY  SOLUTION 
2.1  Method  of  Solution 

A  time-harmonic  e  iwt  plane  wave  E-polarized  in  the  y  direction  is  incident 
broadside  along  the  negative  z-axis  of  a  perfectly  electrically  conducting  (PEC) 
cube  shown  in  Figure  1.  The  front  and  back  faces  of  the  cube  are  located  at  z  =  0 
and  z  =  -s,  the  top  and  bottom  faces  at  y  =  +s  and  y  =  -s,  and  the  right  and  left  side 
faces  at  x  =  +s  and  x  =  -s,  respectively.  To  obtain  a  high  frequency  solution,  we 
begin  by  dividing  the  currents  excited  on  the  surface  of  the  cube  into  the  physical 
optics  (PO)  current  on  the  front  face  and  the  remaining  or  "fringe  currents" 
emanating  from  the  edges  of  the  cube.  Separation  of  the  surface  currents  into  PO 
currents  and  fringe  currents  conforms  to  the  approach  of  Ufimtsev.  5  Although 
fringe  currents  generally  become  negligible  within  a  small  fraction  of  a  wavelength 
from  an  edge,  a  notable  exception  occurs  for  grazing  H-wave  incidence  like  that 
across  the  top  and  bottom  faces  of  the  cube  in  Figure  1,  where  the  fringe  current 
does  not  approach  zero  with  increasing  distance  from  the  leading  edge. 


=  H0»x« 
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Figure  1.  Perfectly 
Conducting  Cube  of 
Side  s  With  Broadside 
Plane-Wave  Incidence 


5.  Ufimtsev,  P,  la.  (1957)  Approximate  computation  of  the  diffraction  of  plane  electro¬ 
magnetic  waves  at  certain  metal  bodies,  Sov.  Phys.  -  Tech.  Phys. ,  27i8):  1708-17  18. 


7 


The  PO  surface  current  is  defined  as  2  e  x  H.  ,  that  is,  the  current  on  the 

z  me 

front  face  of  the  cube  that  would  be  excited  by  the  incident  field  if  the  front  face 
were  part  of  an  infinite  PEC  plane. 

For  the  determination  of  singly  diffracted  rays  from  the  front  edges,  the  fringe 

currents  along  the  leading  front  four  edges  of  the  cube  are  approximated  by  the 

fringe  currents  (total  current  minus  PO  current)  that  would  emanate  from  each  edge 

if  it  were  part  of  the  corresponding  infinite  right-angled  wedge  illuminated  by  the 

incident  plane  wave.  In  the  context  of  GTD,  a  similar  technique  for  evaluating  edge 

diffracted  fields  in  the  vicinity  of  a  caustic  region  has  been  called  the  method  of 

6  7 

equivalent  (edge)  currents.  ’ 

To  determine  the  singly  diffracted  rays  from  the  back  edges,  the  fringe  current 
emanating  from  the  back  four  edges  of  the  cube  can  also  be  approximated  by  those 
of  the  corresponding  infinite  right-angled  wedges,  provide  the  fields  impinging 
upon  these  edges  can  be  estimated.  The  fields  impinging  upon  the  trailing  edges  of 
the  cube  consist  of  the  fields  scattered  by  the  leading  edges  in  addition  to  the  inci¬ 
dent  plane-wave  field.  The  usual  method  of  estimating  the  high  frequency  fields 
impinging  upon  a  trailing  edge  at  grazing  incidence  uses  the  fields  propagated  along 
the  grazing  face  away  from  the  leading  edge,  under  the  assumption  that  the  leading 
edge  were  part  of  an  infinite  wedge.  For  a  right-angled  wedge  the  grazing  currents 

g 

or  fields  have  a  simple  closed-form  expression  which  approaches  the  incident 
plane-wave  field  for  large  s.  Unfortunately,  for  a  cube  side  length  s  less  than  IX, 
where  meaningful  comparison  between  the  high  frequency  solution  and  an  accurate 
integral  equation  solution  is  possible,  the  front  face  of  the  cube  is  a  poor  approxi¬ 
mation  to  the  infinite  face  of  the  right-angled  wedge  and  the  simple  right-angled 
wedge  expression  for  grazing  fields  becomes  a  poor  approximation  to  the  actual 
surface  fields.  Moreover,  this  basic  high  frequency  approximation  for  broadside 
scattering  from  the  cube  remains  inadequate  regardless  of  how  many  interactions 
(multiply  diffracted  rays)  between  the  edges  of  the  cube  are  included  in  the  solution. 
In  order  to  overcome  this  difficulty  and  enhance  the  high  frequency  solution  to  hold 
with  respectable  accuracy  even  in  the  resonance  region,  we  use  a  magnetic  field 
equation  to  integrate  the  PO  and  fringe  currents  on  the  leading  face  and  edges  to  get 
a  better  approximation  to  the  fields  impinging  upon  the  trailing  edges  of  the  cube. 

6.  Ryan,  C.E.,  Jr.,  and  Peters,  L.  ,  Jr.  ( 1969)  Evaluation  of  edge-diffracted 

fields  including  equivalent  currents  for  the  caustic  regions,  IEEE  Trans. 
Antennas  Propag. ,  A.P- 17(3  ):292  -299. 

7.  Knott,  E.  F.  ,  and  Senior,  T.  B.  A.  ( 1974)  Comparison  of  three  high-frequency 

diffraction  techniques,  Proc.  IEEE,  JKh  1 1):  1468-  1474. 

8.  Bowman,  J.  J.  ,  Senior,  T.B.A.,  and  Uslenghi,  P.L.E.  11969)  Electromagnetic 

and  Acoustic  Scattering  by  Simple  Shapes,  North-Holland,  Amsterdam, 
Chapter  6. 


After  the  backscattered  far  fields  are  obtained  in  this  manner  from  the  PO 
current  and  the  singly  diffracted  rays,  the  high  frequency  solution  is  further  refined 
to  include  the  doubly  diffracted  rays. 

Direct  influence  of  any  fringe  currents  existing  along  the  four  edges  of  the  cube 

parallel  to  the  z-axis  is  neglected.  Diffraction  from  the  corners  of  the  cube,  which 

is  higher  order  than  edge  diffraction,  is  approximated  implicitly  when  the  fringe 
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currents  are  integrated  over  the  finite  length  of  the  edges. 


2.2  The  Physical  Optics  Field 

We  begin  the  enhanced  high  frequency  diffraction  (EHFD)  solution  by  finding 
the  backscattered  far  fields  of  the  physical  optics  surface  current 


KPO  =  2  e 


x  TT.  = 
inc 


2  H  S 

°  y 


(1) 


on  the  front  face  of  the  cube,  where  Hq  is  the  amplitude  of  the  incident  H-field  at 
z  =  0.  Insertion  of  the  PO  current  (1)  into  the  integral  for  far-zone  magnetic  field 
produces,  upon  integration,  the  PO  backscattered  magnetic  far  field. 
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where  k  =  2ir/A  is  the  propagation  constant.  Any  valid  far-field  solution  to  the  cube 
should  approach  the  PO  result  [Eq.  (2)],  the  PO  backscattered  magnetic  far  field. 


2.3  The  Singly  Diffracted  Fields  From  the  Front  Edges 

To  obtain  the  singly  diffracted  far  fields  of  the  fringe  currents  of  the  four  front 

5 

edges  of  the  cube,  we  use  the  Ufimtsev  diffraction  coefficients  f  and  g  (for  a  right- 
angled  wedge)  multiplied  by  the  conversion  factor 


which  accounts  for  the  finite  length  s  of  the  edges  of  the  cube.  (This  factor  is 
merely  the  ratio  of  the  far  field  of  a  line  source  of  length  s  to  that  of  an  infinite 
line  source.  ) 


9.  Sikta,  F.A.,  et  al  (1983)  First-order  equivalent  current  and  corner  diffraction 
scattering  from  flat  plate  structures,  IEEE  Trans.  Antennas  Propag. 

.AP-3 1(4):584  -589. 

10.  Collin,  R.  E. ,  and  Zucker,  F.  J.  (1969)  Antenna  Theory,  Part  1,  McGraw-Hill 
New  York,  Eq.  (2.23a). 
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Specifically  for  the  plane-wave  incidence  shown  in  Figure  1,  the  H  far  field 
produced  in  the  backscatter  direction  by  the  fringe  currents  of  the  edges  at 
(x  =  ±s/2,  z  =  0)  is  found  from  the  E-wave  coefficient  f  to  be 


H 

vyfi  ° 


^ikz 

kz 


(4a) 


The  subscript  1  denotes  singly  diffracted  fields  and  the  superscripts  x  and  o  denote 
the  x  edges  on  the  front  face  z  =  0. 

Similarly,  the  singly  diffracted  backscattered  H-field  from  the  fringe  currents 
of  the  front  edges  of  the  cube  at  (y  =  ±  s/2,  z  =  0)  is  found  from  the  H-wave  coeffi¬ 
cient  g: 


H  y°  =  ks 
1  7,3^3 
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aikz 

kz 


e  . 
x 


(4b) 


2.4  The  Singly  Diffracted  Fields  From  the  Back  Edges 


The  far  fields  scattered  by  the  back  edges  of  the  cube  can  also  be  found  simply 
from  the  diffraction  coefficients,  provided  the  grazing  field  impinging  upon  the  back 
edges  of  the  cube  is  known.  Consider  first  the  top  and  bottom  rear  edges  at 
(y  =  ±s/2,  z  =  -s).  Initially  assume  an  average  magnetic  field  impinging  upon  these 
edges  given  by  exp(iks),  but  postpone  the  evaluation  of  until  the  later  Sec¬ 
tion  2.  6.  A  convenient  way  to  determine  the  far  fields  of  these  back  two  edges  is  to 
note  that  only  the  fringe  current  on  the  back  faces  of  these  edges  contributes  to  the 
backscattered  far  field;  thus  the  desired  backscattered  H  far  field  will  be  just  the 
negative  of  the  forward  scattered  field.  Since  this  forward  scattered  H-field  is 
easily  found  from  the  H-wave  diffraction  coefficient  g,  the  backscattered  H  far  field 
from  the  top  and  bottom  back  edges  is  readily  determined  as 


ks  IT,  „2iks  eikz  ^ 

- —  H  e  —t —  e 

7t3\^3  ° 


(5) 


2iks 

The  e  phase  factor,  present  in  (Eq.  (5)]  but  missing  in  (Eqs.  (4)],  accounts  for 
the  extra  propagation  path  length  to  and  from  the  back  of  the  cube.  The  diffracted 
field  [Eq.  (5)]  can  also  be  found  from  the  more  conventional  GTD  procedure  of 
taking  half  the  value  of  the  GTD  H-wave  diffraction  coefficient  when  the  direction 
of  propagation  lies  along  a  terminated  conducting  plane.  11  To  the  field  in  [Eq.  (5)] 
could  be  added  the  secondary  scattered  field  of  the  current  excited  on  the  front  face 


11.  Kouyoumjian,  R.G.  (1975)  The  geometrical  theory  of  diffraction  and  its  applies 
tion,  in  Numerical  and  Asympototic  Techniques  in  Electromagnetics  (Topics 
in  Applied  t’hysics.  Volume  3),  R.  Mittra,  Ed.,  Springer-Verlag,  New  York 
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of  the  leading  edges  by  the  return  field  diffracted  from  the  back  edges.  This  doubly 
diffracted  field  will  be  determined  a\ong  with  the  other  doubly  diffracted  rays,  in 
the  next  section. 

The  grazing  field  travelling  toward  the  left  and  right  rear  edges  of  the  cube  is 
greatly  diminished  by  propagation  along  the  perfectly  conducting  side  faces  parallel 
to  the  electric  field  of  the  incident  plane-wave.  In  other  words,  for  the  plane-wave 
incidence  of  Figure  1,  the  side  faces  severely  shadow  the  back  edges  at 
(x  =  ±s/2,  z  =  -s)  to  the  extent  that  diffraction  from  these  edges  can  be  neglected. 

To  prove  that  this  is  indeed  the  case  even  for  ks  on  the  order  of  1,  we  use  Eq.  (6.  13) 
of  Reference  8  to  estimate  the  current  impinging  upon  the  rear  side  edges,  then 
Eq.  (6.  10)  of  Reference  8  to  find  the  diffracted  H-field  returning  from  an  equivalent 
grazing  plane  wave  that  would  produce  the  impinging  current,  and  finally  the  f- 
diffraction  coefficient  to  find  the  diffracted  far  field  from  the  front  side  edges  subject 
to  the  returning  grazing  H-field.  In  all  we  find  an  H  far  field  diffracted  from  the 
rear  side  edges  given  approximately  by 


/IieiW4 

\  7T  &I<i  H0 


2iks  e 


~TF  ex 


Because  is  less  by  over  an  order  of  magnitude  than  the  far  field  TT^S  diffracted 
from  the  top  and  bottom  back  edges  when  ks  is  greater  than  unity,  and  because  it 
becomes  increasingly  smaller  with  increasing  ks  than  even  the  doubly  diffracted 
terms  (derived  in  the  next  section),  H^s  will  not  be  included  in  the  EHFD  solution. 


2.5  The  Doubly  Diffracted  Fields 

The  doubly  diffracted  rays  are  incident  rays  that  diffract  from  two  edges  before 
returning  to  the  backscatter  direction.  For  the  cube  under  broadside  incidence, 
they  consist  of  three  sets  of  rays:  (1)  those  that  travel  once  between  opposite  edges 
of  the  front  face,  (2)  those  that  travel  once  between  opposite  edges  of  the  back  face, 
and  (3)  those  secondary  rays  mentioned  in  Section  2.  4  that  are  excited  at  the  front 
edges  by  the  rays  returning  from  the  corresponding  trailing  edges.  To  obtain  the 
field  of  each  doubly  diffracted  ray,  the  applicable  diffraction  coefficient  is  used 
twice— the  first  time  to  determine  the  diffracted  field  radiating  from  the  first  edge 
toward  the  second  edge,  and  the  second  time  to  evaluate  the  backscattered  field 
radiated  by  the  second  edge  under  the  illumination  by  the  fields  from  the  first  edge. 
Since  the  three  sets  of  doubly  diffracted  fields  are  evaluated  by  a  straightforward 
manipulation  of  the  diffraction  coefficients,  only  the  results  will  be  given  here. 

Letting  H^0,  TI^S,  anc*  ^ \ Z  designate  the  H-wave  magnetic  field  of  the  three 
sets  of  doubly  diffracted  rays  numbered  above  as  (1),  (2),  and  (3),  respectively, 
we  find 
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The  double-diffraction  factor  A  in  Eqs.  7a,  7b,  7c  equals  2/(37;  ksteiI7^. 
Similarly,  designate  the  doubly  diffracted  E-wave  magnetic  fields  byH^0, 
anc*  The  anc*  from  the  doubly  diffracted  rays  can  be  elimi¬ 

nated  immediately  because  the  singly  diffracted  fields  exciting  these  rays  propagate 
across  the  faces  x  =  ±s/2  parallel  to  the  E-field  and  were  shown  in  Section  2.4  to 
have  negligible  magnitude  for  ks  greater  than  about  1.  The  field  TTg0  of  the  E-wave 
rays  diffracted  across  the  front  face  of  the  cube  from  the  edge  at  x  =  +  s/2  to  the 
edge  ±s/2  has  the  considerably  higher  value. 


tjX°  _  A 
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However,  TT*0  a*so  not  be  included  in  the  EHFD  solution  because  it  decays  with 
ks  faster  than  the  H-wave  doubly  diffracted  terms  [Eq.  (7)],  and  computations  show 
that  it  has  negligible  effect  on  the  final  solution  even  for  ks  smaller  than  1. 


2.6  Evaluation  of  the  Fields  Impinging  Upon  the  Top 
and  Bottom  Trailing  Edges 

The  scatterered  far  fields  given  by  Eqs.  (5),  (7b),  and  (7c)  require  the  total 
magnetic  field  exp(iks)  illuminating  the  top  and  bottom  trailing  edges  of  the  cube. 
As  a  first  approximation  to  one  might  use  the  value  of  the  total  field  leaving  the 
top  and  bottom  leading  edges  of  the  cube  assuming  that  each  of  these  leading  edges  is 
part  of  an  infinite  right-angled  wedge.  Specifically,  the  solution  to  the  right-angled 

g 

wedge  reveals  that  this  total  grazing  field  behaves  asymptotically  (with  large 
distance  ks  from  the  edge'  as 

Ho  ~  Ho(1+  <9) 


For  a  right-angled  wedge,  the  expression  of  Eq.  (9'  for  grazing  field  is  reasonably 
accurate  even  for  ks  on  the  order  of  unity.  However,  its  useful  accuracy  for  the 
cube  is  limited  to  much  larger  values  of  ks.  This  is  because  the  approximation  in 
Eq.  (9)  to  Hq  becomes  accurate  for  the  cube  only  when  the  PO  surface  current  on 
the  front  face  of  the  cube  extends  over  a  large  enough  area  to  produce  the  grazing 
fields  of  the  infinite  front  face  of  the  right-angled  wedge. 
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To  get  a  more  accurate  approximation  for  Hq  than  that  of  Eq.  (9),  we  use  an 
integral  formulation  for  the  magnetic  field  in  terms  of  surface  current  on  the  cube. 
In  particular,  the  total  magnetic  field  [Hq  exp(iks)ex]  at  a  distance  ks  along  the 
center  of  the  top  or  bottom  surface  (x  =  0,  y  =  ±  s/2,  z  =  -s)  of  the  cube  for  the 

broadside  plane-wave  incidence  of  Figure  1  is  expressed  conveniently  by  the  mag- 
3 

netic  field  equation  as 
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where  K  is  the  surface  current,  ,  r)  is  the  scalar  Green's  function 

exp(ik|  r '  -  "?)/!’?'  -  7| ,  and  Sq  is  the  surface  of  the  cube  excluding  the  top  or 

bottom  face  on  which  7q  is  located.  Because  Hq  refers  to  the  magnetic  field 

impinging  upon  the  top  or  bottom  trailing  edge  of  the  cube,  the  current  FT  includes 

only  the  PO  and  fringe  currents  of  the  front  face  and  leading  edges.  In  other  words, 

exp(iks)  is  the  total  magnetic  field  that  would  exist  at  the  surface  point  7q  of  the 

cube,  if  the  sides  of  the  cube  continued  to  infinity  past  the  trailing  edges. 

The  major  contribution  from  the  integral  in  Eq.  (10)  is  from  the  PO  current  on 

PO 

the  front  face  of  the  cube.  Denoting  this  integral  byT  ,  that  is. 
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we  find  upon  substitution  of  iTPO  from  Eq.  (1)  into  Eq.  (11)  that 

iksa^{<!>) 

iks  ,  .  ,  2  -iks  “  i 
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The  function  £*o(0)  is  defined  as 


^  1  +  secd(<t,/4) 
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The  double  integration  of  Eq.  (11)  reduces  to  the  single  integration  of  Eq.  (12a) 
after  converting  the  x',  y^  integration  variables  to  spherical  coordinates 

(R1  sin  9  cos  q,  R1  sin  9  sin  p)  centered  at  the  observation  point  r  (so  that  the 

2  '  0 
area  element  equals  R'  sin  6  d0  d0  and  R  equals  s/cos  9).  and  performing  the  6 

integration.  The  remaining  p  integration  of  Eq.  (12a)  is  done  numerically  in  a  short 

subroutine  of  the  computer  program  for  evaluating  the  EHFD  solution. 

PO 

In  addition  to  the  PO  current  K  ,  the  fringe  currents  associated  with  the  four 
leading  edges  of  the  cube  contribute  to  the  integral  in  Eq.  (10)  for  the  field  impinging 
upon  the  top  and  bottom  trailing  edges.  However,  a  rough  estimate  of  the  contribution 
from  the  top  or  bottom  front  edge  on  the  opposite  back  edge  shows  that  it  is  small  com  - 
pared  with  the  integral  in  Eq.  (12a),  and  thus  can  be  neglected  in  the  determination  H^. 

To  obtain  the  magnetic  field  impinging  upon  the  top  or  bottom  trailing  edge 
from  the  fringe  current  of  its  respective  leading  edge,  we  note  that  for  an  infinite 
right-angled  wedge  <*•**>,  the  integral  in  Eq.  (12a)  becomes  negligible,  and  the 

H  contribution  of  the  PO  current  reduces  to  simply  H  .  Thus  the  A/3  term  in  the 

o  t  r  J  o 

large  ks  expression  in  Eq.  (9)  for  Hq  must  be  produced  by  the  fringe  current  of  the 
respective  leading  edge.  Labelling  this  contribution  to  the  integral  in  Eq.  (10)  by 
Tv,  we  have 


The  contribution  of  Eq.  (13)  to  the  field  impinging  upon  a  trailing  edge  from  the 
leading -edge  fringe  current  can  also  be  obtained  from  the  fringe  diffraction  coeffi¬ 
cient  g  by  realizing  that  the  magnetic  field  radiated  by  the  fringe  current  toward 
the  trailing  edge  is  minus  twice  the  magnetic  field  radiated  in  the  opposite  direction. 

The  final  contribution  to  the  current  integral  in  Eq.  (10)  is  from  the  fringe 
currents  along  the  front  side  edges  at  (x  =  ±s/2,  z  =  0).  At  first  thought  one  might 
argue  that  these  side  fringe  currents  contribute  negligibly  to  the  fields  on  the  top 
and  bottom  faces  because  they  radiate  like  cylindrical  waves  over  the  side  faces. 
Indeed  this  is  true  for  large  ks,  but  for  s  <  A/2  (ks  it)  the  fringe  currents  on  these 
side  edges  will  also  radiate  over  the  top  and  bottom  faces  of  the  cube.  We  can 
estimate  the  contribution  of  the  front  side  fringe  currents  to  the  integral  in  Eq.  (10) 
by  judiciously  applying  the  Ufimtsev  diffraction  coefficient  f.  Because  the  field  of 
the  fringe  current  of  the  corresponding  right-angled  wedge  is  analytic  except  in  the 
direction  of  the  face  of  grazing  incidence,  the  diffraction  coefficient  f  is  analytic 
throughout  the  same  region  and  will  give  the  field  produced  by  the  fringe  current 
inside  the  90°  wedge  angle  as  well  as  outside.  Use  of  this  property  of  f  allows  the 
integral  in  Eq.  (10)  of  the  fringe  currents  along  the  front  side  edges  to  be  approxi¬ 
mated  by 
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H  e 
o  x 


where  f  ,  the  E-wave  diffraction  coefficient  in  the  direction  of  r  ,  equals  approxx- 
o  u 

matety  1/2.  The  spreading  factor  Fg  is  given  approximately  by 


•/i  (-‘i)'-'"4  ( 


sin  (ks/4) 


and  accounts  for  the  finite  length  and  angular  displacement  of  the  radiating  edges. 
The  sin(ks/4)/(ks/4)  part  of  the  spreading  factor  occurs  twice— once  to  account  for 
the  field  point  rQ  on  the  top  or  bottom  face  lying  at  an  angle  roughly  equal  to 
0.  5  radians  off  boresight  of  the  radiating  edges,  and  a  second  time  to  account  for 
the  linear  phase  variation  of  the  radiated  field  across  the  trailing  edges.  (The  first 
portion  of  the  spreading  factor,  that  d  +  )  e  ^  in  Eq.  (14b),  is  an 

asymptotic  result  that  holds  for  ks  <  2;  however,  because7x  has  negligible  effect 
on  the  total  solution  for  ks  >  2,  regardless  of  whether  the  exact  or  asymptotic  form 
of  Eq.  (14b)  is  used,  the  simpler  asymptotic  result  in  Eq.  (14b)  is  retained  for  all 
values  of  ks.  ) 

In  all,  is  found  by  adding  the  individual  current  contributions  of  Eqs.  (12), 
(13),  and  (14)  that  comprise  the  integral  in  Eq.  (10): 

jr/2  iksa0W 

H'=HF'=H  l+-e~ikS  /  - - - —  dt  +  F  A  +  £  .  (15) 

o  o  o  o  V  a_  (0)  s  ’  3 


1  +  l  e'ikS 
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The  factors  A,  (0),  and  Fg  are  defined  in  Eqs.  (7),  (12b),  and  (14b),  respectively. 

Note  that  for  large  ks.  the  amplitude  Hq  of  the  magnetic  field  impinging  upon  the 

top  and  bottom  trailing  edges  reduces  to  Hq,  the  asymptotic  amplitude  of  Eq.  (9)  of 

the  grazing  magnetic  field  propagating  over  the  top  of  the  corresponding  right-angled 

wedee,  Usine  H  '  instead  of  H  in  our  high  frequency  solution  distinguishes  the 

”o  °  2411. 

EHFD  solution 'from  an  ordinary  HFD  solution  or  the  usual  GTD  solution,  *  ’  and, 

as  the  next  section  shows,  extends  the  useful  range  of  application  of  the  high  fre¬ 
quency  solution  to  frequencies  in  the  resonance  region. 
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2.7  Numerical  Results  for  the  High  Frequency  Solutions 
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Now  that  the  physical  optics  far  field  has  been  determined  in  Section  2.  2,  the 
singly  diffracted  far  fields  in  Sections  2.  3  and  2.  4,  the  doubly  diffracted  far  fields 
in  Section  2.5,  and  the  grazing  field  Hq  illuminating  the  top  and  bottom  trailing  edges 
in  Section  2.6,  the  EHFD  magnetic  far  field  (TT^*36'  backscattered  by  the  cube  under 
the  broadside  plane- wave  incidence  of  Figure  1  can  be  written 


TTcube  =  e  H 

SC 


ks 
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(16) 


The  physical  optics,  singly  diffracted,  and  doubly  diffracted  fields  are  given  re¬ 
spectively  by  the  first  second,  and  third  term  within  the  brackets  of  Eq.  (16).  The 
"double-diffraction  factor"  A  is  defined  under  Eqs.  (7),  and  the  "enhancement 
factor"  F^  is  defined  by  Eq.  (15). 

The  radar  cross  section  0 ,  which  is  related  to  by  the  simple  expression 


4*z2  |TT, 


a  = 


cube  I  2 
sc 


(17) 


H 


is  computed  using  the  EHFD  solution  of  Eq.  (16)  and  plotted  in  Figure  2  as  a  func¬ 
tion  of  the  perimeter  (4s /A)  of  the  cube.  Comparison  of  Figure  2  with  the  method 
of  moments  (MOM)  solution  for  the  magnetic  field  integral  equations  of  Section  3 
(Figure  7)  shows  that  the  two  RCS  curves  agree  closely  for  1.5  <  4s/A  <  6,  and 
even  reasonably  well  in  the  resonance  region  (4s /A  <  1.  5)  where  a  high  frequency 
diffraction  solution  is  not  expected  to  yield  valid  results. 

Figure  3  reveals  that  this  unexpected  agreement  in  the  resonance  region  occurs 
because  is  used  instead  of  Hq  for  the  fields  impinging  upon  the  top  and  bottom 
trailing  edges.  The  curve  in  Figure  3  is  the  RCS  computed  from  the  less  accurate 
high  frequency  approximation  of  using  Hq  for  the  grazing  field  at  the  trailing  edges. 
The  variations  of  this  conventional  HFD  solution  in  Figure  3  are  less  pronounced 
then  the  EHFD  solution,  the  resonance  is  hardly  discemablev  and  the  distinction 
between  the  conventional  HFD  solution  and  the  PO  solution  shown  by  the  curve  in 
Figure  4  is  considerably  less  than  between  the  EHFD  and  PO  solution. 
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Figure  2.  Broadside  RCS  Computed  From  Enhanced  High 
Frequency  Diffraction  Solution  (EHFD) 
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Figure  3.  Broadside  RCS  Computed  From  Conventional 
High  Frequency  Diffraction  Solution  (HFD) 
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Figure  4.  Broadside  RCS  Computed  From  Physical 
Optics  (PO)  Solution 
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Figure  5  plots  only  the  PO  plus  single-diffraction  terms  of  the  EHFD  solution 
[Eq.  (16)]  in  order  to  show  the  importance  of  the  contribution  of  the  doubly  diffracted 
rays.  Expressions  for  the  far  fields  of  the  triply  diffracted  rays  were  also  deter¬ 
mined,  but  are  not  shown  because  they  had  no  perceptible  effect  on  either  the  EHFD 
curve  of  Figure  2  or  the  conventional  solution  curve  of  P'igure  3.  This  negligible 
effect  of  the  triply  diffracted  rays  on  broadside  RCS  supports  the  conclusion  that  the 
conventional  high  frequency  solution  is  limited  by  the  use  of  Hq  for  the  fields 
grazing  upon  the  trailing  edges  and  not  by  the  neglect  of  multiply  diffracted  rays  of 
order  higher  than  two. 

3.  THE  MAGNETIC-FIELD  INTEGRAL  EQUATION  SOLUTION 

The  MFIE  has  been  solved  numerically  for  the  broadside  RCS  of  the  cube  by 

1  12 
Tsai,  Dudley,  and  Wilton,  and  more  recently  by  Wang  and  Drane,  up  to  a  cube 

sidelength  of  about  0.  7A.  Unfortunately,  the  MFIE  (and,  incidentally,  the  electric- 

field  integral  equation  EFIE)  suffers  from  spurious  homogeneous  solutions  at  the 


12.  Wang,  J.  J.  H. ,  and  Drane,  C.  J.  (1982)  Numerical  analysis  of  arbitrarily 
shaped  bodies  modeled  by  surface  patches,  IEEE  Trans.  Microwave 
Theory  and  Tech. ,  J£TTV30<8):  1167-  1173. 


1_ 


18 


» 


0  I  2  3  4  6  0  7  0  S  10  II  12  13  14  16 


4  s/X 


Figure  5.  Broadside  RCS  Computed  From  EHFD  Solution 
Without  Doubly  Diffracted  Rays 


frequencies  of  the  corresponding  interior  cavity.3,  13,  14  The  first  interior  resonant  ^ 

frequency  of  the  cube  occurs  at  s/A  =  0.7,  and  subsequent  resonant  frequencies  are 
spaced  at  increasingly  narrower  frequency  intervals.  Consequently,  for  s/A  greater 
than  about  0.5,  numerical  solutions  to  the  MFIE  rapidly  deteriorate.  Figure  6  re¬ 
veals  the  severity  of  the  spurious  resonance  problem  by  plotting  the  numerical 

MFIE  solution  for  RCS  out  to  a  cube  side  length  of  about  4.  The  numerical  scheme  ^ 

that  was  used  is  identical  to  that  of  Reference  1  with  384  patches  over  the  cube.  The 
jagged  peaks  and  valleys  in  the  RCS  curve  of  Figure  6,  are  caused  by  the  spurious 
solutions  of  the  MFIE  at  the  interior  resonant  frequencies. 

13.  Murray,  F.  H.  ( 1931)  Conductors  in  an  electromagnetic  field,  (E°e^>t,  H°e^>t), 

Am.  J.  Math.,  5 3( 2) : 27 5- 28 8.  • 

14.  Maue,  A.W.  (1949)  On  the  formulation  of  a  general  scattering  problem  by 
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Figure  6.  Broadside  RCS  Computed  From  MFIE  Solution 
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A  number  of  methods  have  been  devised  to  eliminate  the  spurious  resonances 
for  the  MFIE  (and  EFIE)  solution  (see  Reference  3  for  a  survey  of  the  various 

3 

methods).  A  particularly  simple  technique  for  eliminating  the  spurious  solution 
"augments"  the  original  magnetic-field  integral  equation  with  the  condition  that  the 
normal  H-field  be  zero  at  the  surface  of  the  perfectly  conducting  scatterer.  This 
augmentation  eliminates  the  spurious  resonances  without  sacrificing  the  basic 
simplicity,  solution  capability,  or  pure  magnetic-field  character  of  the  original 
MFIE.  Specifically,  the  augmented  magnetic-field  integral  equation  (AMFIE)  takes 
the  form 


=  k  "  x  K(r) 


K  x  V  \j/ ds'  . 


(18) 


where  K(r)  is  the  surface  current  at  an  observation  point  7  on  the  surface  s  of  the 
cube,  TI.J7)  is  the  incident  magnetic  field  at  that  observation  point,  and  i p  is  the  usual 
scalar  Green's  function  defined  after  Eq.  (10).  The  original  MFIE  is  reclaimed 
from  the  AMFIE  merely  by  crossing  the  normal  h  to  the  surface  at  r  into  both  sides 
of  Eq.  (18).  Like  the  original  MFIE,  the  integral  in  Eq.  (18)  is  evaluated  in  a 

3 

surface  principal  value  sense. 

A  desirable  feature  of  the  AMFIE  (Eq.  (18)]  is  its  solvability  by  the  same  simple 
numerical  "method  of  moments"  scheme  applied  to  the  original  MFIE.  In  particular. 
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when  we  apply  the  numerical  scheme  that  Tsai,  Dudley,  and  Wilton1  used,  that  is, 
the  same  numerical  scheme  that  we  applied  to  the  MF1E  to  obtain  Figure  6,  the 
HCS  curve  shown  in  Figure  7  results.  (Actually,  the  MFIE  is  used  for  4  s/A  <1.5 
and  the  AIV1F1E  for  4  s/A  >1.5  because  the  MFIE  yields  slightly  greater  accuracy 
at  the  very  low  frequencies.  )  One  can  see  from  Figure  7  that  the  spurious 
resonances  of  Figure  6  are  removed.  However,  the  solution  beyond  4  s/X  *  6  is 
not  shown  because  numerical  experimentation  with  patch  size  showed  that  accept¬ 
able  accuracy  could  be  assured  only  up  to  4  s/X  =  6  when  limited  to  a  patch  size 
corresponding  to  384  patches  over  the  cube.  Smaller  patch  sizes  were  not  possible 
with  the  present  computer  program,  because,  even  using  the  symmetry  of  the  cube 
to  reduce  computer  time  and  storage,  the  central  memory  of  the  available 
Cyber  750  was  not  large  enough  to  handle  more  patches.  Also,  more  than  two  hours  of 
computer  time  was  required  to  obtain  the  curve  of  Figure  7.  (Such  large  computer 
time  and  storage  requirements  are  the  major  reasons  for  the  surface  integral 
equations  having  seen  little  application  to  general  three-dimensional  bodies  greater 
than  a  couple  of  wavelengths  across.  ) 


Figure  7.  Broadside  RCS 
Computed  Using  AMFIE  Solution 
Beyond  4  s/X  =  1.5 


4.  THE  COMPOSITE  SOLUTION 


The  chief  advantage  of  the  numerical  integral  equation  solution  is  its  high 
accuracy  at  low  frequencies  from  the  Rayleigh  region  through  the  resonance  region. 
Its  chief  disadvantage  is  the  prohibitive  amount  of  computer  time  and  storage  that 
limits  its  reliable  accuracy  only  to  the  low  frequency  range  plotted  in  Figure  7. 

Conversely,  the  chief  advantage  of  conventional  high  frequency  diffraction 
solutions  is  their  simplicity  of  expression  and  the  inconsequential  amount  of  com¬ 
puter  time  and  storage  they  require  to  compute  RCS  at  an  arbitrarily  high  frequency. 
Their  major  disadvantage  is  their  restriction  to  frequencies  well  above  resonance 
for  useful  engineering  accuracy.  Consequently,  for  three-dimensional  scatterers 
there  usually  exists  an  intermediate  frequency  range  between  resonance  and  sub¬ 
stantially  higher  frequencies  in  which  neither  the  integral  equation  solution  nor  the 
high  frequency  diffraction  solution  yields  accurate  results.  Moreover,  such  an 
intermediate  frequency  range  in  which  neither  solution  is  accurate  prevents  reliable 
comparison  between  the  integral  equation  and  high  frequency  solutions. 

Fortunately,  we  have  been  able  to  enhance  the  conventional  high  frequency  solu¬ 
tion  for  the  cube  in  order  to  reduce  the  limits  of  accuracy  to  a  fraction  of  a  decibel 
at  frequencies  bordering  on  the  dominant  resonance.  In  addition,  use  of  the  aug¬ 
mented  magnetic  field  integral  equation  has  extended  the  accurate  range  of  our 
integral  equation  solution  to  frequencies  further  above  resonance  than  previous  inte¬ 
gral  equation  solutions  for  broadside  RCS  of  the  conducting  cube.  In  short,  the  EHFD 
solution  of  Figure  2  and  the  AMFIE  solution  of  Figure  7  compare  very  closely  over 
the  intermediate  frequency  range  from  4  s/A  *  1.6  to  6,  and  thus  can  be  combined  to 
form  a  composite  solution  that  holds  for  all  frequencies.  Figure  8  plots  this  compo¬ 
site  solution  out  to  4s/X  equal  to  15  and  compares  this  theoretical  curve  with 

measured  RCS  data  taken  at  11  X-band  frequencies  on  18  solid  aluminum  cubes  vary- 
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ing  in  side  length  from  5  to  100  millimeters.  The  composite  curve  in  Figure  8 
uses  the  AMFIE  solution  of  Figure  7  for  4  s/A  <  5.  1  and  the  EHFD  solution  of 
Figure  2  for  4  s/A  >  5. 1.  The  predicted  and  measured  values  of  RCS  compare  to 
within  the  two-sigma  experimental  limits  of  accuracy  of  about  ±  1  dB. 

Finally,  for  the  sake  of  comparison  with  a  familiar  canonical  solution.  Figure  9 
shows  the  RCS  of  a  perfectly  conducting  sphere  of  radius  "a"  computed  from  the 
exact  spherical  eigenfunction  solution  (for  example,  see  Reference  8,  Chapter  9) 
and  plotted  to  the  same  scale  as  Figure  8.  The  first  resonance  of  the  sphere  occurs 
at  nearly  the  same  perimeter  and  has  nearly  the  same  magnitude  as  that  of  the  cube, 
but  for  frequencies  beyond  the  first  resonance  the  strong  specular  reflection  from  the 
front  face  of  the  cube  raises  its  RCS  well  above  the  RCS  of  the  sphere. 

15.  McGahan,  R.V.  (1983)  Scattering  Experiments  at  the  Ipswich  Electromagnetic 
Measurements  Facility:  Backscatter  FVom  Metal  Cubes.  RA bC  In-House 
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Figure  8.  Broadside  RCS  Computed  Using  Composite  of  AMFIE 
and  EHFD  Solutions.  Circles  indicate  measured  data 
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Figure  9.  RCS  for  Perfectly  Conducting  Sphere  of  Radius  a 
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5.  SUMMARY  AND  CONCLUSIONS 

The  sparsity  of  accurate  (-  i  1  dB)  broadband  solutions  to  general,  three- 
dimensional  scatterers  with  edges  and  corners  motivated  our  looking  into  the  prob¬ 
lem  of  broadside  backscattering  from  the  perfectly  conducting  cube.  Existing  mag¬ 
netic-field  integral  equation  (MFIE)  solutions  for  the  cube  encounter  spurious 
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resonances  that  cause  serious  inaccuracies  for  cube  side  lengths  greater  than  about 
0.  5  A.  Existing  GTD  solutions  for  the  cube  yield  an  accurate  broadside  RCS  only 
at  considerably  higher  frequencies  where  the  GTD  solution  practically  coincides 
with  the  physical  optics  integration  from  the  front  face  of  the  cube.  Thus,  not  only 
is  there  an  intermediate  range  of  frequencies  in  which  no  accurate  solution  exists, 
but  there  is  also  no  common  range  of  frequencies  that  can  be  used  to  compare  reli¬ 
ably  the  low  frequency  integral  equation  solutions  and  the  high  frequency  diffraction 
solutions. 

The  present  approach  has  been  first  to  develop  an  enhanced  high  frequency  dif¬ 
fraction  (EHFIJ)  solution  that  determines  the  RCS  of  the  cube  with  considerable  accu¬ 
racy  down  to  frequencies  in  the  resonance  region;  and  secondly,  to  use  an  augmented 
magnetic -field  integral  equation  (AMFIE)  to  eliminate  the  spurious  resonances  and 
extend  the  range  of  validity  of  the  original  MFIE  up  to  frequencies  well  above  resonance, 
that  is,  to  as  high  a  frequency  as  available  computer  time  and  storage  permits.  This 
twofold  approach  of  using  the  AMFIE  and  EHFD  solutions  enabled  us  to  make  a 
mutually  confirming  comparison  of  RCS  computed  by  two  very  different  techniques 
over  an  intermediate  range  of  frequencies.  Moreover,  it  produced  an  accurate  compo¬ 
site  expression  for  the  broadside  RCS  of  the  cube  at  any  frequency  that  compared  to 
within  ±  1  dB  with  measured  data  taken  for  the  cube  side  length  ranging  from  0.  15 
to  4  wavelengths. 

The  general  development  and  numerical  verification  of  the  augmented  integral 
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equations  were  presented  in  previous  publications.  However  the  technique  used 
to  enhance  the  conventional  high  frequency  diffraction  solution  is  apparently  new, 
although  quite  simple  in  essence  and  straightforward  in  application.  Instead  of  using 
the  conventional  GTD  approximation  for  the  grazing  fields  impinging  upon  the  trailing 
edges  of  the  cube,  we  computed  the  scattered  part  of  the  impinging  fields  by  inte¬ 
grating  over  the  physical  optics  and  fringe  currents  of  the  front  face  and  edges  of 
the  cube.  In  other  words,  the  currents  induced  on  the  leading  face  and  edges  of  the 
cube  were  used  to  better  estimate  the  fields  impinging  upon  the  trailing  edges.  This 
refinement  increased  considerably  the  accuracy  of  the  solution  for  cube  side  lengths 
less  than  one  wavelength,  and  produced  a  high  frequency  diffraction  curve  that  agreed 
closely  with  the  integral  equation  solution  at  intermediate  frequencies.  Unlike  the 
conventional  high  frequency  solutions,  the  EHFD  solution  clearly  revealed  the 
dominant  resonance  of  the  cube  that  peaks  at  a  cube  side  length  of  approximately 
0.  25  wavelengths. 

Finally,  this  enhancement  technique  can  also  be  applied  to  perfectly  conducting 
scatterers  other  than  the  cube,  provided  a  definite  separation  exists  between  leading 
and  trailing  ends  of  the  scatterer  along  the  direction  of  propagation  of  the  incident 
field.  Whether  the  technique  is  feasible  for  more  general  geometries  and  aspect 
angles  remains  to  be  determined. 
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